Selectivity in binary fluid mixtures: static and dynamical properties 
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Selectivity of particles in a region of space can be achieved by applying external potentials to 
influence the particles in that region. We investigate static and dynamical properties of size selec- 
tivity in binary fluid mixtures of two particles sizes. We find that by applying an external potential 
that is attractive to both kinds of particles, due to crowding effects, this can lead to one species 
of particles being expelled from that region, whilst the other species is attracted into the region 
where the potential is applied. This selectivity of one species of particle over the other in a localized 
region of space depends on the density and composition of the fluid mixture. Applying an external 
potential that repels both kinds of particles leads to selectivity of the opposite species of particles 
to the selectivity with attractive potentials. We use equilibrium and dynamical density functional 
theory to describe and understand the static and dynamical properties of this striking phenomenon. 
Selectivity by some ion-channels is believed to be due to this effect. 

PACS numbers: 82.70.-y, 05.20.Jj, 61.20.Gy, 87.16.dp 



I. INTRODUCTION 



Biological ion channels are amazing nanofluidic de- 
vices. Ion channels are special proteins with pores 
through the center that allow for passive transport of 
ions, such as K + and Na + , along their electrochemical 
gradients, through a membrane formed by a lipid bilaycr. 
The pores of these proteins have a diameter of order a few 
Angstroms and, in general, have a charge of their own. 
Ion channels have two important functions: (i) They can 
open and close the pore and thereby control the current 
through the channel. This phenomenon is called gating, 
(ii) Typically, they can select the type of ions that can 
pass through the pore, a phenomenon called selectivity. 
Selectivity can occur with respect to several properties. 
Some channels can select divalent ions over monovalent 
ones. Other channels can distinguish ions of the same va- 
lency by size. Using gating and selectivity, ion channels 
are responsible for a wide range of physiological phenom- 
ena such as the regulation of ion concentrations inside the 
cell [1]. 

The selectivity with respect to ions of equal electri- 
cal charge (e.g., between K + and Na + ) is rather puz- 
zling. Recent theoretical studies of ion channels with 
wide pores, such as the 1-type Ca channel, examined the 
effects of entropy (also referred to as molecular crowding) 
in the three component mixture composed of neutral sol- 
vent molecules and two ion species, which differ only in 
size, and also studied the influence of the electrostatic 
attraction of the ions into the channel [2]. It was found 
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that the entropic effects can dominate over the electro- 
static attractions. The theories used in these studies were 
(i) equilibrium density functional theory (DFT) for mix- 
tures of hard spheres in an external field and (ii) bulk 
fluid models of charged hard spheres, in which the ex- 
ternal potentials are mapped on to shifts in the chemical 
potentials of the different components [3, 4, 5]. 

In recent years, colloidal suspensions have become pop- 
ular model systems for testing all aspects of liquid state 
theory, including DFT. In contrast to atomic or molecu- 
lar fluids, colloidal suspensions can conveniently be stud- 
ied optically, e.g., by confocal microscopy. In addition, 
individual or groups of colloids can be manipulated using 
optics: particles which are optically denser than the sur- 
rounding solvent are attracted to regions of higher light 
intensity, e.g, to a laser focus. This is the working prin- 
ciple of optical tweezers, which allows for the creation of 
an attractive potential well for the colloids, by means of 
laser irradiation. Because of this ability to manipulate 
and observe the individual colloids, we believe that, as 
has been the case for other liquid state phenomena, col- 
loidal suspensions provide a useful model system in which 
to understand how selectivity in some ion channels oc- 
curs. Thus, although the present study is motivated by 
the desire to further understand entropy driven selectiv- 
ity in ion channels, we consider a more general scenario. 
Using an accurate DFT, we study equilibrium selectivity 
and using dynamical density functional theory (DDFT) 
we study the dynamics of selectivity in colloidal binary 
mixtures. 

Dynamical aspects of selectivity have neither been 
studied experimentally nor theoretically, up to now. 
While a microscopic dynamical theory for so-called 'sim- 
ple' liquids is still under construction [6, 7, 8, 9, 10], the 
DDFT for Brownian particles (i.e. for colloidal suspen- 
sions) [11, 12, 13] has been applied with success to a large 
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FIG. 1: We consider a binary colloidal mixture confined in a 
three-dimensional slit pore of width 2 L. In the center of the 
slit pore an attractive or repulsive external potential V ex t acts 
on the particles. This external potential can lead to selectivity 
of one species of particle over the other in this region of the 
slit. 



number of situations, ranging from phase separation [14] 
and spinodal decomposition [15] to the micro-rheology of 
colloid polymer mixtures [16, 17]. 

Here, we use DDFT to study the dynamics of selectiv- 
ity in a binary colloidal mixture of particles of two differ- 
ent sizes. We consider the situation where the mixture 
is confined within a three-dimensional slit pore, where, 
apart from the usual fluid ordering near the walls of the 
slit, the equilibrium fluid is homogeneously distributed 
within the slit. At time t = 0, we switch on an exter- 
nal potential in the center of the channel (see Fig. 1) and 
then follow the time evolution of the system using DDFT. 
We select this configuration since it may be realized ex- 
perimentally in a straightforward manner, by using laser 
tweezers to create the external potentials. 

This paper is laid out as follows: In Sec. II we give a 
short overview of static and dynamical DFT and describe 
the model system under consideration. Sec. Ill begins 
with a discussion of the influence of external potentials 
on the equilibrium fluid density distributions. This is 
followed by a discussion of the dynamics of selectivity 
for a binary mixture of hard spheres. We close with a 
summary and an outlook in Sec. IV. 



II. THEORY AND MODEL SYSTEM 

The theories that we use for describing the fluid are 
static (equilibrium) DFT and DDFT. Here, we only in- 
troduce the aspects of the theories that are relevant to the 
present study. For a detailed introduction to DFT, see 
e.g. Rcfs. [18, 19, 20] and for an introduction to DDFT, 
see Refs. [11, 12, 13, 15]. 

A. Equilibrium DFT 

We consider a fluid mixture composed of v different 
types of particles at a temperature T, confined within a 
fixed volume T. For systems where external potentials 



V^t( r ) arc acting on particles of type z, with i = l...u, 
treating the system in the grand canonical ensemble, one 
can rigorously prove the existence of the grand potential 
(free energy) functional f2[{/?j}]. The set of equilibrium 

fluid density profiles {p^ eq ^(r)}, are those which minimize 
f2[{pi}]. The minimum value of f2[{/Jj}] is the grand po- 
tential of the system [20] . The grand potential functional 
may be written in the following way: 



d 3 r 



(1) 



where fii is the chemical potential for particles of type i 
and the intrinsic Hclmholtz free energy functional 

F[{f>i}\ = E \ J t lnA ?^( r ) - !] d 3 r + F ex [{ Pi }}, 

(2) 

which is the sum of the ideal-gas parts, where 1/0 = 
ks T is the thermal energy, Ai is the thermal Dc Broglie 
wavelength, and ./^[{pi}] is the free energy contribution 
originating from the interactions between the particles. 

From the minimization principle on f2[{pi}], it follows 
that the equilibrium fluid density profiles are the solution 
of the following Euler-Lagrange equations: 

mPl}] = = r 1 lnAfp i( r) + ^Ml +y «( r) _ M , 



6pi(r) 



(3) 

For a given set of external potentials {V^t(r)}, these 
yield a set of non-linear equations for the density pro- 
files {pi(r)}, which can be solved e.g. using an iterative 
numerical algorithm. 

The excess free energy functional ^^[{jOj}] is only 
known exactly for a few one-dimensional model systems. 
However, for many experimentally relevant systems, ap- 
proximate functionals have been constructed which yield 
results that are in quantitative agreement with experi- 
ments and simulations. In particular, this is the case 
for mixtures of hard spheres, for which functionals aris- 
ing from fundamental measure theory (FMT), mainly 
based on geometric considerations, have been developed 
[21, 22, 23]. In FMT the excess free energy functional is 
given by 



Pf ex [{pi}] = / $({n Q (r)})d 3 r, 



(4) 



where the excess free energy density $ is a function of a 
set of weighted densities 



a (r)=E / ft(r>i(r- rOdV. 

i=l ^ 



(5) 



In Eq. (5) the geometrical weight functions are denoted 
by u>a(r) and a labels four scalar and two vector- like 
weight functions. It is the fundamental measure theory 
constructed in Ref. [22] that we implement in this study. 
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B. Dynamical DFT 

In a manner analogous to the Runge-Gross theorem for 
quantum systems [24], one can prove the existence of a 
dynamical density functional theory for classical systems 
[25]. As in the static case, the proof of the existence of a 
DDFT does not lead to having a theory with which one 
can do calculations. However, for systems of Brownian 
particles with over-damped stochastic equations of mo- 
tion [see Eq. (6), below], a successful DDFT has been 
developed by Marconi and Tarazona [11, 12] for calcu- 
lating the time evolution of the ensemble average density 
profiles pi(r,t), i = 1 . . . v. Note that by 'ensemble av- 
erage' density, we mean an average over the ensemble 
of different realizations of the stochastic noise terms and 
particle starting positions in the equations of motion. For 
a detailed discussion of the differences between the dy- 
namics of the ensemble averaged, the instantaneous, and 
the coarse grained density, see Ref. [13]. 

In real three-dimensional systems, the number of parti- 
cles in a volume T can be changed only by a flux through 
the boundaries cT, which, in absence of chemical reac- 
tions, directly leads to local particle conservation. In 
Sec. II A above, we discussed DFT, an equilibrium sta- 
tistical mechanical theory for non-uniform fluid mixtures 
in the grand canonical ensemble. However, for systems 
with impermeable boundaries, the equilibrium limit leads 
to the canonical ensemble. For this reason we consider 
fixed numbers Ni of particles of type i at positions r$ (t) , 
where n = 1 . . . Ni, and we assume that the particle dy- 
namics is governed by the following over-damped Brow- 
nian equations of motion: 



dr 



(') 



dl 



7i 



W«(r«) + FWl + v^NWft), 



(6) 

(i) 

where F„ is the net force on particle n of type i, due to 
the interactions with all the other particles in the system 
and 7i is the mobility coefficient for particles of type i. 

is a Gaussian random white noise, with zero mean 
and the correlator 



NW(t)oN&>(f; 



2S nm 6 ij 6(t-t')l. (7) 



The Fokker-Planck equation corresponding to Eq. (6) 
describes the time evolution of W({r^}, t), the probabil- 
ity of finding the particles in the system at the positions 
r^} at time t. By integrating over W({r { n ] },t) with re- 
spect to all but one of the coordinates of particles of type 
i, one obtains a time evolution equation for the ensem- 
ble averaged density pi(r, t) [14]. For systems of particles 
interacting via pairwisc interaction potentials, the result- 
ing time evolution equations for pi (r, t) depend only on 
the non-equilibrium two-body distribution functions, but 
for systems interacting via many-body interaction poten- 
tials, the time evolution equations for pi(r,t) depend on 
the higher-body distribution functions [15]. The time 
evolution equations for pi(r, t) are the lowest members in 



a hierarchy of equations that is similar to the 

Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hier- 
archy [6, 7, 13, 14, 18]. The DDFT is obtained by using a 
closure relation to truncate the hierarchy. This is done by 
approximating the non-equilibrium two-body and higher- 
body distribution functions by the same quantities in an 
equilibrium system with the same one-body density dis- 
tributions [11, 12, 15]. One can show that for a system 
with a given set of density distributions {pi(r)}, one can 

find a unique set of external potentials {V^t( r )} 5 such 
that the system would be at equilibrium if it was ex- 
posed to these external potentials. This closure leads to 
the following set of coupled DDFT equations [14]: 



dpi(r,t) 
dt 



7i V 



Pi(r,t) V 



Spi 



(8) 



where the functional f2[{/?j}] is the equilibrium grand po- 
tential functional, given in Eq. (1). Note, that the gradi- 
ent of the variation of the grand potential equals the gra- 
dient of the variation of the total free energy, i.e., the sum 
of the intrinsic free energy .F[{/5;}] and the contributions 
from the external potentials V ext - For non- interacting 
particles, Eq. (8) reduces to the drift-diffusion equation. 
The (equilibrium) closure approximation used to obtain 
these equations means that there are some limitations 
on what kinds of problems the DDFT can be applied 
to. For example, one can not use it to describe barrier 
crossing in systems with free energy barriers. However, 
some aspects of glassy systems can be described within 
the DDFT framework [26]. 

In order to implement Eq. (8) using the FMT approx- 
imation for f2[{pi}], we may take advantage of certain 
simplifications that arise due to the structure of the FMT 
when calculating the gradient of the variation of the ex- 
cess free energy. We present this in the Appendix. 



C. Model system 

The model system that we study is a binary mixture of 
hard-spheres, composed of small (s) and big (b) particles 
with sphere radii R s and i?b, respectively. The particle 
densities are p s (r, t) and Pb(r, t). We first consider a bulk 
mixture with given densities and examine the influence 
of external potentials, which can be cither attractive or 
repulsive, that act on a localized region. 

Depending on the origin of the external potentials, 
they may act in the same way on both types of parti- 
cles, or alternatively the potentials may be proportional 
to the volume of the particles. In studies of ion channels 
with wide selectivity filters with fixed charges, such as the 
1-type calcium channel, it was found that a competition 
between energy and entropy can explain size selectivity 
[3, 4, 5]. In such ion-channels the external potentials act- 
ing on the ions are generated by the electric field due to 
the charges fixed in the selectivity filter. Since the po- 
tentials generated in this way are the same for particles 



4 



system 1 (reservoir) 




system 2 


ii! 1 ': fixed 




ii ( V-V '".-fixed 

' ! r I CXI 




exchange 




<- 


of particles 




p.'".- fixed 




p ; ; calculated 



FIG. 2: The effect of size selectivity can be rationalized in 
a simple model of two coupled systems. In system 1, which 
acts as a reservoir for system 2, the densities of the mixture 
are fixed to be p^. In system 2, the chemical potentials are 
shifted relative to the corresponding values in system 1 by the 
external potentials V}*\, which are considered to be constant 
throughout the system. The different particle densities in 
system 2 result from these shift in the chemical potentials. 



carrying the same charge, we model these by setting the 
potentials acting on both species of particles to be equal, 



i.e. we set V}^(r) = V^l(r), to model a mixture with 
two sizes of particles having the same sign and magni- 
tude charge on them. 

However, one may also consider a different situation: 
that of e.g. optical tweezers applying a force on col- 
loidal particles within a specified region of space. When 
the forces acting on a colloidal mixture are generated 
by optical tweezers, then the external potentials act- 
ing on the particles are proportional to the volume of 
the particles. To model this situation, one should set 



r{b) 



V c ^(r) = (R s /R b fV^(r). 

The mechanism of the size selectivity as studied in the 
next section can be understood by considering a simpli- 
fied system of a binary colloidal fluid mixture partitioned 
into two parts: in the first part, a finite region of space, in 
which the external potentials arc applied, and the second 
part, which acts as a 'reservoir', consisting of the remain- 
ing parts of the system. If we assume both systems to be 
infinite in size (i.e., we neglect interface and wall effects) 

(i) 

and if we assume V^ xt to be constant in system 2, we ar- 
rive at the situation depicted in Fig. 2 with the chemical 
potentials and p,^ 2 ' in system 1 and 2, respectively 
given by 



0), 



p 



5T_ 

5pt 



(9) 



for i = s, b and £=1,2 and with the intrinsic free energy 
JF[{pi}] defined in Eq. (2). py is the spatially constant 
density of species j in system i. 

In system 1 (which is the reservoir) we specify the den- 
sities of the mixture, p^\ where i = s,b, which fixes the 
chemical potentials in system 1. In system 2, which is 



1.1 



1.05 



0.95 



0.9 



/ 

_J 


\ 


— small spheres ; 

— big spheres > 

■ 


i 

1 



10 



15 



20 

z/R s 



25 



30 



35 



FIG. 3: Density profiles pi(z\ where i = s,b, divided by 
the respective bulk densities p\ ulk (corresponding to the bulk 
packing fraction r\i = |j?f7rpi ), of a binary mixture of 
hard spheres with a size ratio of Ri, = 2R S . The density 
profiles show small-particle selectivity: locally, with the help 
of the external potentials, the density of the small spheres 
(solid line) is slightly increased, while the density of the big 
spheres (dashed line) is slightly decreased. Using Eq. (11), it 
is possible to calculate the external potentials V^ xt (z), i = s,b, 
which give rise to these equilibrium density profiles. These 
external potentials are displayed in Fig. 4, for various bulk 
packing fractions. 



coupled to system 1, so that particle exchange with the 
reservoir is allowed, we apply the (spatially constant) ex- 
ternal potentials V^ t . Applying these external potentials 
is then equivalent to shifts in the chemical potentials in 

(2) 

system 2. The densities in system 2, p\ , can thus be 
computed by solving the following set of equations 



£\{p { i ) })=»T\{pf ) }) + V£ t for i = s,b. 



(10) 



If the external potentials arc attractive, then the chemical 
potentials in system 2 are effectively increased over those 
in system 1, while they are decreased in the case of repul- 
sive potentials. In the low density limit, the Eqs. (10) for 
i = s and i = b decouple into a pair of independent equa- 
tions. At higher densities, Eqs. (10) form a non-linear 
pair of coupled equations for the densities in system 2. 
As discussed in Sec. IV below, these equations allow us to 
understand the phenomena described in the next section. 
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FIG. 4: In order to generate a slight local increase in the den- 
sity of the small spheres and a slight decrease in the density of 
the big spheres (as displayed in Fig. 3), it is necessary to ap- 
ply external potentials V^ xt {z). Four examples are displayed 
here for a binary mixture with size ratio Rt = 2R 3 and equal 
bulk packing fractions 77s = rjb- At low bulk packing fractions 
f]a — Vb = 0.01, we see in (a) that an attractive potential has 
to act on the small spheres and a repulsive potential on the 
big spheres, in order to generate local small-particle selectiv- 
ity. This can be understood from the ideal-gas limit. In (b), 
where rj s — r/b = 0.04, the repulsion acting on the big spheres 
is significantly reduced, as a result of the competition between 
entropy and energy. In (c), where rj a — rjb — 0.09, both exter- 
nal potentials are roughly the same and are attractive and in 
(d), where rj 3 — rjb — 0.15, the attraction required to attract 
the small particles is weaker than the attraction required to 
effectively repel the big particles. 

III. RESULTS 

A. The external potentials 

Within the framework of static DFT, the equilibrium 
thermodynamic properties of the system are obtained 
from the solution of the Euler-Lagrange equations (3), 
obtained by minimizing the functional in Eq. (1). Typ- 
ically, one specifies the external potentials V^t( r ) an d 
then computes the resulting equilibrium density distri- 
butions Pi(r), i = 1, . . . , v. However, it is also possible 
to employ the Euler-Lagrange equations (3) in order to 
compute the external potentials that give rise to a speci- 
fied set of density distributions {pi(r)} [27]. One obtains: 

« = _ ^IM _ r i lnA 3 (r) 

5pi(r) 

It is worth noting that in the limit of vanishing densi- 
ties, the contribution due to particle interactions, which 



depends on the variation of J^aiKPi}], becomes negli- 
gible and one recovers the ideal-gas result: V^t( r ) = 
— /3 _1 In Af jOj(r) + fii. In the ideal-gas limit, where the 
density distributions are simply the Boltzmann factors 
of the external potentials, it is clear that in order to lo- 
cally increase the density of the particles over the bulk 
value, attractive potentials are required, and to locally 
decrease the density, repulsive potentials are needed. 

For a one-component interacting system, the behav- 
ior is similar to the ideal-gas case, although the density 
distribution is not simply the Boltzmann factor of the ex- 
ternal potential. The density can locally be increased or 
decreased by an attractive or repulsive external potential, 
respectively. 

The situation becomes more interesting in the case of 
a binary mixture. If we start with constant densities 
in an infinite bulk system and wish to generate a set of 
density profiles such as those shown in Fig. 3, which show 
a local increase in the density of the small spheres and 
at the same location a decrease in the density of the big 
spheres, we can use Eq. (11) to calculate the external 
potentials required to achieve this. 

In Fig. 4 we display the external potentials that one 
must exert on the system in order to observe the two 
density profiles shown in Fig. 3, for a binary hard-sphere 
mixture with Rt = 2R S and equal bulk packing fractions 
of the small and large particles, rj s and rjb, respectively 
(rji = | i?f 7T Pj Ulk , where p\ ulk is the density of species 
i in the bulk, at points in space a large distance from 
the support of the external potential). At low values of 
Vs = Vb = 0.01, we see in Fig. 4(a), that the required 
external potentials are similar to those predicted by the 
ideal-gas functional. The external potential that gives 
rise to a local increase in the density of the small spheres 
is attractive (solid line) and the external potential that 
gives rise to a decreased in the density of the big spheres 
is repulsive (dashed line). However, if the bulk pack- 
ing fractions increase to r] s = rjb = 0.04, as shown in 
Fig. 4(b), we find that the external potential acting on 
the big particles is significantly less repulsive than in the 
low density case displayed in Fig. 4(a). This stems from 
the competition between energy, due to the interactions 
of the particles with the external potentials, and entropy, 
due to the inter-particle interactions and the resulting ex- 
cluded volume. The influence of the competition between 
energy and entropy is more pronounced if we increase the 
bulk packing fractions even further. For rj s = rjb = 0.09, 
in Fig. 4(c) we see that both species of particles have to 
be subjected to an attractive potential well, in order to 
obtain density profiles that locally have an increase in 
the density of the small particles, at the same time as a 
decrease in the density of the big particles. For this par- 
ticular choice of packing fractions, the external potentials 
for both sizes of particles are roughly the same. In order 
keep the reduction of the density of the large particles 
at the same moderate level at much higher packing frac- 
tions rj s = rib = 0.15, the attraction on the big particles 
has to be much stronger than the attraction on the small 
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FIG. 5: The external potential acting on the small spheres 
(solid lines) and big spheres (dashed lines). In addition to 
the soft slit potentials, we perturb the density distributions 
by additional potentials in the middle of the slit, which can 
be either (a) attractive or (b) repulsive. 
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FIG. 6: (a) The equilibrium density profiles of a binary hard- 
sphere mixture with Rb = 2 R 3 confined within a slit with 
potentials given by Eq. (12). (b) When the attractive poten- 
tials in the middle of the slit are switched on [see Fig. 5(a) and 
Eq. (13)], the density of the small spheres is locally enhanced 
in the middle of the slit while the density of the big spheres is 
decreased. This is an example of small particle selectivity. In 
both (a) and (b) the chemical potentials /j,b and fis are chosen 
such that the adsorptions given by Eq. (14) are TbR 2 = 0.495 
and T S R 2 = 1.308. The chemical potentials in (b) correspond 
to reservoir packing fractions r\ a = rj s = 0.15, as indicated by 
the dotted line. 



particles. 



B. Equilibrium DFT for selectivity 

From the results presented so far, we can see that the 
competition between energy and entropy in mixtures can 
lead to surprising effects. The question that arises is 
whether or not it is possible to use this competition in 
order to select one species over the other with the aid 
of an appropriately chosen external potential. To this 
end, we consider a binary hard-sphere mixture confined 
between two planar walls - i.e. in a slit. We model the 
potentials due to the walls of the slit as follows: 



V&(z) = A 



(12) 




(b) 
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FIG. 7: Same as Fig. 6, but for repulsive potentials, (a) The 
equilibrium density profiles with Rb = 2 R a confined within 
a slit with potentials given by Eq. (12). (b) When the re- 
pulsive potentials in the middle of the slit are switched on 
[see Fig. 5(b) and Eq. (13)], the density of the small spheres 
is locally reduced in the middle of the slit while the density 
of the big spheres is increased. This is an example of large 
particle selectivity. In both (a) and (b) the adsorptions given 
by Eq. (14) are T b R 2 b = 0.786 and T s R 2 a = 0.905. The chem- 
ical potentials in (b) correspond to bulk reservoir fractions 
rj a — r] s — 0.15, as indicated by the dotted line. 



where A = 10 k B T, z = 20.48 R s , m = 20 and L l = 
L + Ri, where L = 15 R s . Due to the large value of the 
power to, Eq. (12) models the continuous potentials of 
a pair of slightly soft parallel planar walls, separated a 
distance 2L. zo is the center of the slit pore. When the 
fluid mixture in the slit is at equilibrium, we perturb the 
fluid by introducing additional external potentials in the 
middle of the slit, so that the external potentials are: 



V&{z) = A 



Z - Z 

u 



e» exp 



z - zo 
w 



(13) 



The additional second term in the external potentials 
can be either attractive (ej < 0), as shown in Fig. 5(a), 
or repulsive (ej > 0), as indicated in Fig. 5(b). The range 
of the external potential is w = 3 R s and the depth or 
height 6j = ±2 kgT. 

First, we consider the case of an attractive potential 
in the middle of the slit (e, < 0). We compare two equi- 
librium states: (i) the binary mixture in the slit with- 
out the additional potentials (ej = 0), and (ii) the equi- 
librium state of the binary mixture with the potentials 
(f-i = —2ksT), shown in Fig. 5(a). We calculate the den- 
sity profiles which minimize the equilibrium DFT by solv- 
ing the Eulcr-Lagrange equations, Eq. (3). Since equilib- 
rium DFT is defined in the grand canonical ensemble, 
particle exchange between the system and a reservoir, 
which sets the chemical potentials, is possible. However, 
we wish to compare our equilibrium DFT calculations 
with DDFT computations, which preserves the number 
of particles. In order to make an honest comparison, we 
minimize the equilibrium DFT under the constraint that 
the amount of material within the slit is fixed, i.e., that 
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the adsorptions 



T t = / Pi (z)dz, (14) 

J — oo 

are fixed. To enforce these constraints, we treat the 
chemical potentials /i Q and /!& as Lagrange multipliers 
and we choose them such that in the case when the po- 
tentials in the middle of the slit are switched on, the reser- 
voir packing fractions are rj s = r/b = 0.15. In the case of 
an attractive potential this results in T S R 2 S = 1.308 and 
T^R 2 = 0.495, while in the repulsive case the correspond- 
ing adsorptions arc T S R 2 = 0.905 and T^R 2 = 0.786. 
Correspondingly, the reservoir packing fractions when 
the potentials are switched off are different. 

When the attractive potential wells are not switched 
on, we find that the density profile of the big spheres 
close to the walls of the slit possess prominent peaks [see 
Fig. 6(a)], indicating strong packing and correlation ef- 
fects at the walls of the slit. The density of the small 
spheres is slightly decreased at the wall. In the middle 
of the slit we find a region where both densities are con- 
stant. When the attractive potentials in the center of the 
slit are switched on, so that both the small and the big 
particles are subjected to the external potentials shown 
in Fig. 5(a), we find that there is a significant local in- 
crease in the density of the small particles in the region 
of the attraction, while the density of the big particles 
is strongly decreased here - see Fig. 6(b). Hence, if the 
bulk densities of the big and the small particles are suf- 
ficiently large, an attractive potential well can be small 
particle selective due to the competition between energy 
and entropy. 

We now turn our attention to the case when a repul- 
sive potential is turned on in the middle of the slit - see 
Fig. 5(b). In Fig. 7(a) we display the equilibrium density 
profiles for the case when the repulsive potentials are not 
switched on. The profiles are similar to those shown in 
Fig. 6(a). The difference is due to our treatment of the 
chemical potentials as Lagrange multipliers to set the ad- 
sorptions in Eq. (14), for the case when the potentials in 
the middle of the slit are switched on. This leads to 
the adsorptions and therefore also to the density profiles 
being different in the two cases when the potentials are 
switched off. 

In Fig. 7(b) we display the equilibrium fluid density 
profiles for the case when the repulsive potentials in the 
middle of the slit, for both components, are switched 
on. In the region where the repulsive external potentials 
are applied, the local density of the small spheres is de- 
creased and the density of the big spheres is increased. 
It is interesting to note that despite the repulsion acting 
on the big particles, the local density of the big spheres 
is greater than the density in the surrounding fluid. This 
is an example of big-particle selectivity. 



C. Dynamics of selectivity 

In order to better understand the behavior of the sys- 
tem, we study the time evolution of the density profiles 
of both types of particles as they evolve between the two 
equilibrium states shown in Fig. 6, the case with the at- 
tractive potentials applied in the center of the slit, and 
those shown in Fig. 7, the case where the repulsive poten- 
tials are exerted in the center of the slit. In what follows, 
all times t* = tjr are given in units of the Brownian time- 
scale r = ((3/j)R 2 ; it is roughly the time it takes for a 
particle to diffuse over a distance equal to its own radius. 
Note also that we have set the two mobility coefficients 
to be equal, 71 = 72 = 7- In a colloidal suspension this 
would not be the case: the mobility of spherical particles 
in a suspension decreases with the inverse of the cross- 
sectional area. For the systems under consideration here 
this results in ~f s /jb = 4. With the small particles being 
considerably faster than the big ones, the initial behavior 
of the big particles described in the following might not 
be observable in colloidal suspension. 

The initial conditions for attractive and repulsive po- 
tentials in the center of the slit pore are those displayed 
in Figs. 6(a) and 7(a), respectively. 

In the first case, when we switch on the attractive po- 
tentials in the center of the slit at time t* = 0, initially 
both types of particles behave as one would intuitively 
expect: they follow the attraction and move towards the 
center of the slit. However, this is only the case for short 
times. This initial influx of particles results in a slight 
increase in the particle densities at the center of the slit, 
as can be seen in the density profiles shown in Fig. 8. 
For times t* > 1, the drift diffusion behavior of the sys- 
tem qualitatively changes, as the competition between 
the energy due to the interactions of the spheres with 
the external potentials and the entropy due to the hard- 
sphere interactions between particles, sets in. For times 
t* > 1, only the small particles follow the attraction and 
the number density of the small particles in the center 
of the slit increases further, while the big spheres are ex- 
pelled from this region. At short times, only spheres in 
the center of the slit show a net movement and the den- 
sity profiles do not change over time in the vicinity of the 
slit walls. However, at longer times, there is a net flux of 
the small particles from the walls of the slit towards the 
center and there is also a net flux of big particles from the 
center towards the walls. This diffusion process is slow 
and the time it takes for the system to reach the final 
equilibrium state depicted in Fig. 6(b) is rather long. It 
also depends on the system size - the wider the slit, the 
longer it takes for the particles to diffuse (say) from the 
wall to the center. 

We observe similar behavior in the case when we apply 
repulsive potentials at the center of the slit. When we 
switch on the repulsive potentials at t* = 0, both types 
of particles behave as our intuition suggests, and they 
follow the repulsion and move way from the center of the 
slit. This initial flow results in a slight decrease in the 
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FIG. 8: Time evolution of the density profiles of the small (solid lines) and big (dashed lines) spheres from the initial state in 
Fig. 6(a) towards the final state in Fig. 6(b). At time t* = we switch on the attractive potentials in the middle of the slit. 
Note that for short times, both species of particles follow the attraction and move towards the center. However, for larger times 
t* > 1, only the density of the small spheres increases in the region of the potential well, while the big spheres move away from 
the center, despite the attractive external potential. 
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FIG. 9: Time evolution of the density profiles of the small (solid lines) and big (dashed lines) spheres from the initial state in 
Fig. 7(a) towards the final state in in Fig. 7(b). At time t* — we switch on the repulsive potentials in the center of the slit. 
Note that for short times, both species of particles follow the repulsion and move away from the center. However, for larger 
times t* > 1, only the density of the small spheres decreases in central region, while the density of the big spheres increases in 
the center, despite the repulsive external potential. 
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particle densities at the center of the slit, as can be seen in 
the density profiles shown in Fig. 9. However, for longer 
times t* > 1, the behavior of the system qualitatively 
changes, as the competition between the energy and the 
entropy sets in. For times t* > 1, only the small particles 
follow the repulsion and the number density of the small 
particles in the center of the slit decreases further, while 
the big spheres become effectively attracted towards this 
region, and the density of the big particles increases at 
the center of the slit. The final equilibrium state the one 
shown in Fig. 7(b). 

Note that these counterintuitive results occur only 
when the densities of the particles are high enough. At 
low densities, when the system can be modelled as an 
ideal-gas, attractive potentials always lead to an increase 
in the densities of both types of particles and repulsive 
potentials always lead to a decrease. 

IV. DISCUSSION 

Using the bulk approach depicted in Fig. 2 and ex- 
pressed in Eqs. (10) it is possible to understand the re- 
sults presented in Sec. III. If the external potentials v} xt 
are the same for all components of the mixture, as as- 
sumed in our DFT and DDFT calculations, we can re- 
produce the increase or decrease of the particle densities 
in the region where the external potentials are applied: 
for a binary mixture with size ratio Rb = 2 R s and bulk 
packing fractions of ?/ s = rjb = 0.15 in the reservoir sys- 
tem 1 (corresponding to ps = R% = 0.036) and 
an attractive external potential of magnitude 2 ks T, the 
density of the small particles in system 2 increases by a 
factor of 1.91, while the density of the big particles is 
reduced to 0.008 of its original value. This is in good 
agreement with the results displayed in Fig. 6: the den- 
sity of the big particles at the center of the slit pore is 
negligibly small while the density of the small particles 
at the same spot is about twice the reservoir density. In 
the case of a repulsive potential of magnitude 2 ks T, the 
density of the big particles increases by a factor of 2.25, 
while the density of the small spheres is reduced to 0.21 
of its value without the field. Again, this agrees rather 
well with the results shown in Fig. 7. 

While we have shown results for size selectivity only 
for one size ratio, i.e. Rb = 2R S , we have confirmed that 
a binary mixture with less asymmetric radii shows analo- 
gous behavior. However, in the case of a less asymmetric 
size ratio, the amplitude of the attractive or repulsive 
potential has to be larger to generate a degree of size 
selectivity similar to that reported in our study. 

In addition, Eqs. (10) allows one to easily estimate the 
influence of the external potentials in other scenarios and 
determine at what densities the system should or should 
not display selectivity. For example, if the external po- 
tentials are proportional to the volume of the particles, 
which would be the case in an experimental realization 
using laser tweezers, then one should expect results sim- 



ilar to those reported here, for a binary mixture with 
a size ratio Rb = 2 R s and packing fractions r/ s = 0.15 
and r]b = 0.02. In this case, external fields on the order 
of V£l = ±lk B T and v£] = ±8k B T are sufficient to 
generate selectivity. 

In the present study, the relatively wide slit geometry 
has only a small effect on the selectivity observed. How- 
ever, for more narrow pores, the grand canonical calcula- 
tions in a cylindrical channel geometry in Rcf. [2] showed 
that the confinement can enhance the selectivity. 

The arguments presented above are based purely on 
the equilibrium free energy and hold for any type of un- 
derlying dynamics, i.e., for the overdamped Brownian dy- 
namics considered here as well as for systems following 
Hamiltonian dynamics. A key feature of the dynamics 
of separation presented in Figs. 8 and 9 is that the big 
particles, although they are finally driven out of the re- 
gion in which the attractive potential acts and into the 
region in which the repulsive potential acts, initially fol- 
low the direction given by the gradient of the external 
potential. For Figs. 8 and 9 we assumed equal mobilities 
for the big and the small particles. In a colloidal sus- 
pension this is not be the case (the mobility of spherical 
particles in a suspension decreases with the inverse of the 
cross-sectional area) and so the initial behaviour of the 
big particles described above might not be observable in 
a colloidal suspension. 

To conclude, we should remind the reader that the 
DDFT in Eq. (8) that we have used to describe the sys- 
tem, does not incorporate hydrodynamic interactions be- 
tween the colloidal particles. Whilst hydrodynamic inter- 
actions do not affect the static properties of the system 
(i.e. the equilibrium DFT is still applicable), hydrody- 
namic interactions do have an influence on the dynamical 
properties of the system. The influence of hydrodynamic 
interactions between the particles has been incorporated 
in the DDFT in several ways [28, 29]. We believe that ex- 
tending the present study to include the influence of hy- 
drodynamic interactions would not qualitatively change 
any of the results that we observe since these only influ- 
ence the dynamics but not the energetics of the system. 
It is the energetics which determine the final equilibrium 
state, i.e. whether selectivity can be observed or not. 



APPENDIX: DDFT AND THE STRUCTURE OF 
FMT 

For the calculation of the time evolution of the density 
profiles, the gradient of 5fl/5pi(r) are required in Eq. (8). 
This includes the gradient of the variation of the excess 
free energy functional T ex . Since we employ FMT for the 
excess free energy, we can make use of the the structure of 
T exi given in Eqs. (4) and (5). Note that in the general, 
three-dimensional case, the variation of the excess free 
energy with respect to the density profile of component 
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i can be written as 

SJ^ex 



E 

a 

E 



dn a 
dn a 



w i a (i > -r)dV. (A.l) 



It follows from Eq. (A.l) that the gradient acts solely on 
the weight functions: 



V f — 



V^(r'-r) d 



V. 



(A.2) 



Although the sum in Eqs. (A.l) and (A.2) is over four 
scalar (a = 3, ...,0) and two vector-like (a = v2, vl) 
weight functions, one can exploit the relations between 
the scalar weight functions AttR^ uiq = AnRi cu\ = lo\ 
and the vector- weighted functions AirRi w,^ = to l v2 . This 
allows one to reduce to the sum in Eqs. (A.l) and (A.2) 
to three terms by introducing the auxiliary functions 



*«(r) 



<9$ 

dn 3 ' 



1 <9$ 



1 5$ 



dn 2 47ri?i drii 47ri?f dno ' 



and 



*$(r) 



<9$ 



1 <9$ 



dn v2 AirRi dn v i 



(A.3) 



(A.4) 



(A.5) 



These functions depend on the particular version of FMT 
employed. Here we use the White-Bear version [22]. 

The slit geometry that we consider in the present study 
allows us to further simplify the gradient in Eq. (A.2). 
The effective one-dimensional weight functions and their 
derivatives with respect to z are easily calculated. For 
the volume weight function wc find 



ui(z) = ir(R*-z 2 ) Q(Ri-\z\ 



(A.6) 



where O(z) denotes the Heaviside step-function. The 
derivative of the volume weight function is 

Wl(z) = £-u\(z) = -27TZ Q(Ri - |*|). (A.7) 



Note that the derivative of the Heaviside function can be 
neglected, because the weight function <jj\{z) vanishes at 
the integration boundaries. The effective surface weight 
function is 



4(z) = 27ri? 4 0(ifc - \z\), 



(A. 



which is constant in the range of integration, so that its 
derivative 

Wl(z) = ^<4i») = 2irR t [S(z + Ri) - 5(z - R t )\ 

(A.9) 

has contributions only from the derivative of the 
Heaviside-0 function. Finally, the vector-like weight 
function is 

w* 2 (z) - -2irz G(Ri - \z\) e z , (A.10) 
and leads to a derivative 

WUz) = ^«4(*) = I-27T Q(Ri - \z\) (A.ll) 
+ 2TrR q [S(z + R4) + S(z - Ri)]} e z . 

Putting all the ingredients together we obtain the fol- 
lowing expression for the gradient of the variation of the 
excess free energy in planar geometry: 



d 5F ex 
dz Spi(z) 



^\z')Wl{z'-z) 



+ ¥: ] 2 ( Z ') w: 2 ( Z ' - z)} dz' 

+ 2TtR l ^ 2 t \z + R l )-^ ( 2 l \z-R l ) 
+ ^liz + RJ + ^liz-R,)}. (A.12) 

ACKNOWLEDGMENTS 

A. J. A. acknowledges financial support from the British 
Council, funded under the ARC programme, and 
also from RCUK. M.R. acknowledges financial support 
from the priority program SPP 1164 of the Deutsche 
Forschungsgcmcinschaft. M.R. and R.R. acknowledge fi- 
nancial support from DAAD funded under the PPP pro- 
gramme. 



[1] B. Hille, Ion Channels of Exciteable Membranes (Sinauer 

Asc, Sunderland, 2001). 
[2] R. Roth and D. Gillespie, Phys. Rev. Lett. 95, 247801 

(2005). 

[3] W. Nonner, L. Catacuzzeno, and B. Eisenberg, Biophys. 
J. 79, 1976 (2000). 

[4] W. Nonner, D. Gillespie, D. Henderson, and B. Eisen- 
berg, J. Phys. Chem. B 105, 6427 (2001). 



[5] D. Gillespie, W. Nonner, D. Henderson, and R. S. Eisen- 
berg, Phys. Chem. Chem. Phys. 4, 4763 (2002). 

[6] A. J. Archer, J. Phys.: Condens. Matter 18, 5617 (2006). 

[7] A. J. Archer, J. Chem. Phys. 130, 014509 (2009). 

[8] J. G. Anero and P. Espafiol, Europhys. Lett. 78, 50005 
(2007). 

[9] S. Melchionna and U. M. B. Marconi, Europhys. Lett. 
81, 34001 (2008). 



11 



[10] U. M. B. Marconi and S. Melchionna, preprint: 

arXiv:0902.3694 (2009). 
[11] U. M. B. Marconi and P. Tarazona, J. Chem. Phys. 110, 

8032 (1999). 

[12] U. M. B. Marconi and P. Tarazona, J. Phys.: Condens. 

Matter 12, A413 (2000). 
[13] A. J. Archer and M. Rauscher, J. Phys. A: Math. Gen. 

37, 9325 (2004). 
[14] A. J. Archer, J. Phys.: Condens. Matter 17, 1405 (2005). 
[15] A. J. Archer and R. Evans, J. Chem. Phys. 121, 4246 

(2004). 

[16] M. Rauscher, M. Kriiger, A. Dominguez, and F. Penna, 
J. Chem. Phys. 127, 244906 (2007). 

[17] C. Gutsche, F. Kremer, M. Kriiger, M. Rauscher, R. Wee- 
ber, and J. Harting, J. Chem. Phys. 129, 084902 (2008). 

[18] J.-P. Hansen and I. R. McDonald, Theory of simple liq- 
uids (Academic Press, 2006), 3rd ed. 

[19] R. Evans, in Fundamentals of Inhomogeneous Fluids, 
edited by D. Henderson (Marcel Dekker, New York, 



1992), chap. 3, pp. 85-173. 

[20] R. Evans, Adv. Phys. 28, 143 (1979). 

[21] Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989). 

[22] R. Roth, R. Evans, A. Lang, and G. Kahl, J. Phys.: Con- 
dens. Matter 14, 12063 (2002). 

[23] H. Hansen-Goos and R. Roth, J. Phys.: Condens. Matter 
18, 8413 (2006). 

[24] E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 
(1984). 

[25] G. K.-L. Chan and R. Finken, Phys. Rev. Lett. 94, 
183001 (2005). 

[26] A. J. Archer, P. Hopkins, and M. Schmidt, Phys. Rev. E 

75, 040501 (R) (2007). 
[27] J. R. Henderson, Mol. Phys. 74, 1125 (1991). 
[28] C. P. Royall, J. Dzubiella, M. Schmidt, and A. van 

Blaaderen, Phys. Rev. Lett. 98, 188304 (2007). 
[29] M. Rex and H. Lowen, Eur. Phys. J. E 28, 139 (2009). 



